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Abstract 

The baryonic potential in the framework of the SU{'i) random vortex world- 

t:J- , surface model is evaluated for a variety of static color source geometries. For 

^^ \ comparison, carefully taking into consideration the string tension anisotropy en- 

P^ ■ gendered by the hypercubic lattice description, also the A and Y law predictions 

^ , for the baryonic potential are given. Only the Y law predictions are consistent 

^ ' with the baryonic potentials measured. 

>—i ' 

00 
(N 

CN ■ PACS: 12.38.Aw, 12.38.Mh, 12.40.-y 

JT^ ■ Keywords: Center vortices, infrared effective theory, confinement, baryon potential 

(N 

o 
o 
o 

-)— » 
> 

•i-H 

X 



^email: engel@nmsu.edu 



1 Introduction 

The random vortex world- surface model of the Yang-Mills vacuum ^J |21 El E] is designed 
to describe the central phenomena induced by the strong interaction in the infrared 
regime: Confinement, spontaneous breaking of chiral symmetry, and the axial Ua{1) 
anomaly. The motivation for the model is drawn from lattice Yang-Mills studies [3 El 
HI IHl ISl Uni HH 1121 HH mi which suggest the relevance of center vortices for infrared 
strong interaction physics. For the case of SU{2) color, the model has been shown 
to correctly predict (in the sense of agreement with corresponding lattice Yang-Mills 
results) a number of observables characterizing the aforementioned phenomena. 

In a recent initial study of the generalization to SU{3) color ^, the random vortex 
world-surface model was again found to reproduce the confinement properties of Yang- 
Mills theory. Both a low-temperature confining phase as well as a high-temperature 
deconfined phase are generated, separated by a weakly first order phase transition. 
The confinement properties are intimately related to the percolation properties of the 
vortices. 

Apart from the order of the deconfining transition, a further aspect in which the 
SU{3) theory differs substantially from the SU{2) case is the long-range baryonic po- 
tential, which was not studied in [3]. This constitutes the subject of the present work. To 
arrive at definite statements concerning the baryonic potential, two issues must be taken 
into account. One issue is the breaking of spatial rotational symmetry in the present 
hypercubic lattice formulation of the model; after all, to measure baryonic potentials, 
it is necessary to evaluate Wilson loops which are not located within two-dimensional 
lattice planes. On the other hand, the Wilson loop areas spanned by baryonic configura- 
tions are substantially larger than the ones encountered in simple mesonic Wilson loop 
calculations; as a result, it is necessary to employ appropriate numerical noise reduction 
techniques to extract a meaningful signal for the baryonic potential. Having dealt with 
these two issues, the data collected in this work unequivocally point to a Y law for the 
baryonic potential as opposed to a A law, cf. Fig. ^ It should be noted that the form 
of the baryonic potential generated by vortex models has been the subject of recent 
debate [13 dj; on the other hand, also within full SU{3) lattice Yang-Mills theory, 
only quite recently sufficiently large separations have been probed to clearly detect a Y 
law f71dlEll2ni- In the full theory, short-distance perturbative A behavior presum- 
ably masks the long-range Y law such that a crossover between the two dependences 
is observed at intermediate distances JH] and the Y law emerges only at rather large 
separations. Aside from the vortex model, also other models of the strong interaction 
vacuum yield a Y law for the baryonic potential, e.g. the dual superconductor picture 
fn\ 122] and the stochastic vacuum [23] ■ 

2 Model description 

The random vortex world-surface model is based on the notion that the Yang-Mills 
vacuum is populated by random closed lines of quantized chromomagnetic flux (vor- 





Figure 1: Based on the notion that three SU{3) color sources in the Yang-MiUs vacuum are bound 
by (hnear) potentials associated with distances characterizing the spatial geometry of the sources, two 
distinct model bond topologies can be envisaged: The A law, where sources are connected pairwise 
(left), and the Y law, where bonds emanating from the sources meet at an (optimally placed) central 
junction (right). While the string tension in the Y law must be the same as in a mesonic configuration 
(where also a single bond emanates from each color source), the A model is assumed to be associated 
with half that string tension. This is based both on the short-distance perturbative limit and on the 
limit where two color sources approach each other, which should yield the usual mesonic string tension. 
Since in this limit, two of the bonds making up the A configuration merge, each of them is taken to 
contribute half of the mesonic string tension (this linear superposition of bonds constitutes an auxiliary 
model assumption which is not readily justified within the underlying nonabelian Yang-Mills theory, in 
which chromoelectric flux tubes interact nontrivially when approaching each other) . 



tices), described by an ensemble of closed random world-surfaces in four-dimensional 
(Euclidean) space-time. In the present implementation of the model, these random sur- 
faces are generated on a hypercubic lattice, composed of elementary squares on that 
lattice. The action governing the ensemble is related to the surface curvature: If two 
elementary squares which are part of a vortex surface share a lattice link but do not lie 
in the same plane, this costs an action increment c. Formally, this can be written as a 
sum over lattice links, 
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where q^uix) describes the chro mo magnetic flux associated with the elementary square 
extending from the site x into the positive n and u directions. In the SU{3) model, the 
variables q^u{x) can take three values, qfj.u{x) G { — 1, 0, 1}; the value indicates that the 
elementary square is not part of a vortex surface, whereas the values ±1 indicate that it 
is. The latter two values distinguish the two possible types of quantized flux carried by 
SU{3) vortices. These will be characterized in detail in the next section by their effect 
on Wilson loops circumscribing them. Note that S is symmetric with respect to the two 



possible types of vortex flux. Note furthermore that, formally, there are two variables 
connected with the elementary square extending from the site x into the positive // and z/ 
directions, namely q^u{x) and qufiix); these correspond to the two opposite orientations, 
i.e. senses of curl, one may ascribe to a surface. However, since both variables are meant 
to describe the same flux (which, apart from its location in space-time, also carries a 
definite orientation), they are related by^ qtiv^x) = —qu^i{x). 

It should also be emphasized that the vortex ensemble is generated subject to the 
constraint of continuity of flux (modulo 27r, i.e. modulo Dirac strings |21 El El)- For 
one, this forces vortex surfaces to be closed, as already mentioned above; on the other 
hand, the existence of two types of quantized flux in the SU{3) model nevertheless 
allows for vortex branchings, where a vortex of one type splits into two vortices of 
the other type while maintaining continuity of fiux'^. As described in detail in j3], 
continuity of flux in practice is guaranteed during the generation of the vortex world- 
surface ensemble by performing updates simultaneously on the six squares making up 
the surface of an elementary three-dimensional cube in the lattice. This is done in a 
way which corresponds to superimposing the (continuous) flux of a vortex of the shape 
of the elementary cube surface onto the flux previously present. 

The flux world-surfaces defined in this manner are formally infinitely thin; however, 
they are meant to describe physical vortex structures present in the Yang-Mills vac- 
uum which possess a finite transverse thickness, akin to the structures making up the 
Copenhagen vacuum ^^. The model surfaces represent the geometrical centers of such 
physical thick vortices. The thickness enters the model description through the spacing 
of the lattice on which the vortex surfaces are generated. For instance, two parallel 
physical vortices can only be meaningfully distinguished from one another when they 
are further apart than a minimal distance. In the model, this distance is encoded in 
the lattice spacing; at shorter distances, the thick structures cannot be distinguished 
and thus should be represented instead by a single model world-surface of combined 
flux. For a thorough discussion of the role of the lattice spacing in the random vortex 
world-surface model, cf. [Hill. 

Given this definite physical meaning of the lattice spacing, it follows that it is a 
fixed finite quantity in the model^, which should be contrasted with conventional lattice 
gauge theory, where one is ultimately interested in the limit of vanishing lattice spacing. 
This has consequences as far as the breaking of spatial rotational symmetry by the 
underlying hypercubic lattice structure is concerned. Whereas in lattice gauge theory, 
spatial rotational symmetry is restored in the course of taking the continuum limit 

^Note that, in light of the description above, this (correctly) implies that a vortex carrying one type 
of flux is equivalent to an oppositely oriented vortex carrying the other type of flux for the present 
purposes (namely, for evaluating Wilson loops). Indeed, these two cases only differ by a flux 27r, i.e. a 
Dirac string, which is unobservable in Wilson loops. On the other hand, when considering topological 
properties j^ U, the two cases do have to be distinguished and one cannot trade off geometrical 
orientation of a vortex in space-time against the magnitude of flux it carries. 

•^This does not happen in the SU{2) model, which only has one type of flux, corresponding to the 
existence of only one nontrivial element in the center of the SU{2) group. 

"^This should also be expected from the point of view that the random vortex world-surface model 
is an infrared effective model, which is only deflned up to a fixed ultraviolet cutoff. 



PBj . in the random vortex world- surface model the breaking of rotational symmetry 
engendered by the hypercubic lattice description and manifest in the action (Q) is a real 
model property which will be quantified in detail in section Roughly, string tensions 
along diagonal directions in the lattice are enhanced compared with the one measured 
along lattice axes to the extent that the most extreme string tension values obtained are 
separated from their mean by 13 %. This variation must be taken into account when 
comparing baryonic potentials with A or Y law expectations. 

It should be noted that, while the random vortex world-surface model is intended to 
be more than just a qualitative model of strong interaction physics, on a quantitative 
level it is generally not expected to furnish a description of medium- to long-range ob- 
servables better than to within of the order of 10 %. After all, it encodes a truncated 
dynamics in which phenomena such as the Coulomb potential between color sources are 
missing and effects at distances shorter than about 0.4 fm cannot be resolved. In light of 
this, the consequences of the hypercubic realization of the random vortex surfaces, which 
include the breaking of spatial rotational symmetry discussed above, represent just one 
particular model distortion of reality among others. No more (or less) significance needs 
to be attached to this specific distortion than to others on a similar quantitative level, 
such as ambiguities in measurements of the topological susceptibility 2j and of the 
quenched chiral condensate [3]. Note also that the breaking of rotational symmetry is 
only introduced by the particular (technically convenient) hypercubic lattice implemen- 
tation of the model; a manifestly symmetric random vortex world-surface model could 
e.g. be formulated using randomly triangulated surfaces in continuous four- dimensional 
space-time. 

3 Observables 

Physically, the random vortex surfaces represent quantized chromomagnetic flux. This 
means that they contribute in a characteristic way to Wilson loops; if one chooses an area 
spanning a given Wilson loop^, then for each time a vortex world-surface pierces^ that 
area, the Wilson loop acquires a phase factor corresponding to the center of the gauge 
group. For the case of SU{3) color, there are two possible phase factors, associated with 
the two possible types of quantized vortex flux, namely exp(±i27r/3). To be precise, 
consider an elementary Wilson loop (plaquette) extending from y into the positive k 
and A directions, with the integration oriented such that one starts at y, integrates 
first into the positive k, direction, and then onwards around the plaquette. Denote this 
by Uf^xiy). The plaquette Uf^x{y) is pierced precisely by the dual lattice elementary 
square q^u{x), where the indices k, A,/i, z/ span all four space-time dimensions and x = 
1/ + (ck + ca — e^ — ei,)a/2, with a denoting the lattice spacing. U^xiy) can thus be given 

^The choice of area is immaterial due to the continuity of flux. 

^Note that Wilson loops are defined on a lattice dual to the one on which the vortices are defined, 
i.e. on a lattice shifted by the vector {a/2, a/2, a/2, a/2), where a denotes the lattice spacing. Thus, 
the notion of a vortex piercing a Wilson loop area is unambiguous. 



exclusively in terms of the corresponding q^^ (x) , namely 

U^xiy) = exp (m/3 ■ e^Xf,u qf^uix)) (2) 

(with the usual Euclidean summation convention over Greek indices). 

To evaluate an arbitrary Wilson loop, it is sufficient |1| to find a tiling of that Wilson 
loop by a set of plaquettes and multiply the values those plaquettes take in the given 
vortex configuration, as specified by (J2I). 

Consider now a baryonic Wilson loop, 

W = -eabc^a'b'd^T ^2 ^t (3) 

o 

where Fj denotes the path-ordered exponential (Wilson) line integral along the path 
taken by the i-th color source. All three Fj start at the same initial space-time point 
and end at the same final space-time point. Due to the special properties of center 
vortex configurations, the baryonic Wilson loop in such configurations can be straight- 
forwardly decomposed into three standard mesonic Wilson loops. On the one hand, any 
center vortex configuration can be written in terms of gauge fields defined purely within 
the (Abehan) Cartan subgroup of the underlying gauge group |211 . Under these circum- 
stances, the Fj in (jHI) are diagonal, and by multiplying (jH]) with the (unit) determinant 
of another color matrix Fq given by an additional exponential line integral connecting 
the end point of the Fj with their initial point, 

1 = T^^TfVf (4) 

one can rewrite the baryonic Wilson loop as 

^ = \T. 6L(rrrr)(r^'T^^)(F-F-) (5) 

a,h,c 

On the other hand, since in a center vortex configuration, the exponential integral along 
a path with coinciding initial and final points is proportional to the unit matrix (this is 
the defining property of center fiux), this indeed is equal to a product of three standard 
mesonic Wilson loops, 

W^ = Q Tr FiFo) (^ Tr F^Fo) (^ Tr F3F0) (6) 

In practice, in the present work, the baryonic correlator of three Polyakov loops was 
evaluated. This is of course nothing but a special case of the baryonic Wilson loop 
discussed above, with paths F, starting from an initial spatial position xq at time t = 0, 
running to spatial positions Xi while at t = 0, remaining there as time runs through the 
complete extension of the lattice to the final time t = (3 (which is identified with t = 0, 
with periodic boundary conditions on the gauge fields), and returning to xq along the 
same path as they took at t = 0. Due to the Abelian nature of the vortex gauge fields, 
all the line integrals at t = and t = j3 cancel and one is left with a baryonic Polyakov 



loop correlator. In accordance with the discussion above, this is furthermore equal to 
the expectation value of a product of three standard mesonic Polyakov loop pairs such 
as the ones used to evaluate standard mesonic Polyakov loop correlators. The baryonic 
Polyakov loop correlator observable has the advantage that boundary effects resulting 
from the propagation of the color sources from xq to Xj and back cancel; on the other 
hand, it limits the temporal extension of the lattice universe one can achieve without 
losing the baryonic Wilson loop signal to numerical noise^. 

4 Numerical noise reduction 

To nevertheless reduce the numerical noise contaminating the mesonic and baryonic 
Polyakov loop correlator measurements as far as possible, the noise reduction technique 
introduced by Liischer and Weisz J2ZI was employed. Given an action which can be 
decomposed as 

S[q]^S[q\q^] = Sy] + J2Sy,q!], (7) 

i 

an observable O which correspondingly factorizes as 

0[q]^0[q\q']=l[Oy,qt], (8) 

i 

and a constraint 5 which separates as 

S[q]^S[q\q^]=S'[q']ll6y,qt], (9) 

i 

the method of Liischer and Weisz corresponds to evaluating the expectation value of O 
over the set of variables q subject to the constraint 6 using the action S as 

{0)hs;S} = \ \\Dq'\ b\q'\ exp(-S*[g*]) J] /[^g-l'^l?*, ^AO\q\ rf^ exp(-S^[g*, ^^) 

= (]lm{,!;S';S'M) (10) 

\ i I {q;5;S} 

In other words, the variables g* are kept fixed while the inner expectation values are 
taken over their respective variables qf using the actions S^ subject to the constraints 
5*; in the outer averaging over the full set of variables q subject to the full constraint 5 
using the full action S, the product of these inner expectation values may be much easier 
to sample than the original 0[g*, g*] was. Note that, in practice, this method often can 
be iterated in the sense that one may again select from the set of variables g* a subset 
which is to be kept fixed while its complement is being used for averaging; as long as the 
action, the constraint and the quantity being averaged can be decomposed in analogy 

''Corresponding to this treatment of the baryonic potential, also the mesonic string tension will be 
evaluated in section [fusing Polyakov loop correlators. 



to (ITj)-®, one can construct additional hierarchies of averaging, where the quantities 
being averaged at each level only depend on a subset of the variables the quantities at 
the next-lower level depended on. 

This method applies to the random vortex world-surface model if one chooses as the 
variables g* all elementary squares which extend into the time and one space direction, 
and as the variables g| the elementary squares at the i-th lattice time (which extend into 
two spatial directions). The action (^ is a sum over terms associated with lattice links, 
each term coupling only elementary squares attached to the corresponding link. The sum 
over spatial links at the i-th lattice time thus yields the S^ piece of the action, whereas 
the sum over all temporal links yields the S"* piece. Similarly, the constraint of continuity 
of flux is to be satisfied independently at each lattice link ji] by the elementary squares 
attached to that link; it thus factorizes into terms coupling only the variables g* and qf 
at fixed i, as well as a term constraining only the variables g*. Lastly, any Wilson loop 
can be evaluated as a product over single plaquettes, which in turn each only depend on 
the value of a single vortex elementary square according to ^; this can thus be trivially 
grouped into factors satisfying the factorization (jH} . In accordance with the remarks on 
the update procedure in section El in this first level of averaging, elementary updates 
are thus only performed on surfaces of elementary lattice cubes which lie in a fixed time 
slice, i.e., which extend into three spatial directions. 

After this decomposition, in practice also a second averaging hierarchy was employed; 
namely, the set of elementary squares g* extending into the time and one space direction 
was decomposed into sets of squares q^^ connecting the (2z — l)-th lattice time with the 
{2i)-th lattice time, with the remaining elementary squares making up the set g* to be 
kept fixed at this second level of averaging. With this decomposition, the properties 
([Tj)-® are again satisfied; note that, for fixed i, the variables g2j only enter the product 
of inner expectation values 

Thus, in this second level of averaging, updates were performed on the surfaces of all 
elementary lattice cubes except for the ones connecting even lattice times 2z with the 
next higher odd lattice times 2i + 1. 

In practice, it turned out to be efficient to carry out the innermost averaging using 
either 4000 or 8000 configurations; in the second-level averaging, either 200 or 400 
configurations were used. For the outermost averaging, typically 20 to 60 configurations 
were enough to achieve a sufficient level of accuracy. 

5 Angular dependence of the string tension 

As already mentioned further above, the mesonic string tension was evaluated using 
Polyakov loop correlators, employing the noise reduction techniques discussed in the 
previous section. In practice, 16^ x 4 lattices were used, with the curvature coupling c 
in (^ set to the physical value c = 0.21, cf. 0. The temporal extension of four lattice 



spacings corresponds to a temperature of T = 0.45Tc (where Tc denotes the deconfine- 
ment temperature). In the case of the static potential between sources separated along 
one of the lattice axes studied in |1], the string tension at this temperature deviated 
from the extrapolation to zero temperature by less than 1 %. 

The mesonic geometries studied in the present work are listed in Table Q] The ge- 
ometry labels defined there will be used below to refer to the different cases. String 



Geometry 
label 


Relative spatial separation 
of Polyakov loops 


Scale factors 
n used 


String tension in 
lattice units aa^ 


MlOO 


n- (a, 0,0) 


1,2,3,4,5 


0.7659 ±0.0004 


MHO 


n ■ (a,a, 0) 


1,2,3 


0.907 ±0.001 


M210 


n ■ (2a, a, 0) 


1,2 


0.885+^^^ 


M310 


n ■ (3a, a, 0) 


1,2 


0.83+nf 


Mill 


n ■ (a, a, a) 


1,2 


1 n9+tJ.0U2 
-L-U^-0.035 


M211 


n ■ (2a, a, a) 


1,2 


omt'oT 



Table 1: Mesonic geometries and associated string tensions. 



tensions a were extracted by fitting the behavior of the measured Polyakov loop cor- 
relators P{nda) at different spatial distances nda separating the Polyakov loops by the 

ansatz 

— In P{nda) r 



a+- 
n 
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Aanda 

using the ra = 1,2 data, where a denotes the lattice spacing^. Of course, 4a is the 
temporal extension of the lattice, cf. above; thus, the denominator on the left hand side 
is simply the minimal area spanned by the pair of Polyakov loops. The accuracy of 
this procedure was checked by also performing fits of the form a ± r/n ± p/n^ to the 
n = 1,2,3 data in the MlOO and MHO cases; this led to adjustments of the string 
tension values by at most 3 % (in both cases, downwards). These fits then did also 
describe the n = 4, 5 data in the MlOO case within their statistical errors. In the error 
analysis of quantities derived from string tensions extracted using only n = 1,2 data, 
the aforementioned 3 % systematical downwards uncertainty in these string tensions is 
incorporated. The string tensions are listed in Table ^ 

Whereas A law predictions for the baryonic configurations studied further below 
can be made purely on the basis of o"Afioo and cta/iio, to arrive at Y law predictions, 
it is necessary to interpolate continuous angular dependences of the mesonic string 
tension from the discrete string tension data collected. To be specific, it is neces- 
sary to combine ctmioo ) o"m3io ? <^M2io and aMuo to interpolate the string tension for 
two sources separated along the direction (cosa,sina, 0), and to combine <Jmwo , o"M2ii 



^In physical units, obtained 0] by equating the zero-temperature string tension (along lattice axes) 
with (440 Me V)^, the lattice spacing in the random vortex world-surface model is 0.39 fm. 





Figure 2: Angular dependence of the string tension a. Left panel: As a function of the angle a 
when the static sources are separated along the direction (cos a, sin a, 0); right panel: As a function of 
the angle (3 when the static sources are separated along the direction ((l/\/2) sin/3, (1/V2) sin/3, cos/3). 
Curves are fitted to the measured string tension values and to the extreme values allowed by their error 
bars. In the case of cr(/3), the fit to the measured values and the upper extreme coincide within the 
resolution of the plot. 

and (T A/111 to interpolate the string tension for two sources separated along the di- 
rection ((1/v^) sin/3, (1/v^) sin/3, cos/5). Since the functional form of the a- and (3- 
dependences is unknown, the quite conservative corridors for these dependences de- 
picted in Fig. 121 were allowed for. Error bars for Y law predictions below are based on 
the extremes allowed by these corridors. 

6 Baryonic geometries 

Baryonic potentials were measured by evaluating baryonic Polyakov loop correlators as 
discussed in section 01 As in the mesonic case, IQ^ x 4 lattices were used at the physical 
point c = 0.21. Table 121 lists the baryonic geometries studied^. 

Similar to the mesonic case, the leading long-range behavior of the potential was 
extracted by fitting the measured baryonic Polyakov loop correlators P^ at the scale 
factors n by the ansatz 

- In P^ U 

-P±J^ = V + - (12) 

Aana n 

using the n = 1,2 data. The baryonic tension V can thus be directly compared to 

A and Y law predictions obtained by adding bond lengths (in units of no) weighted 

by the associated string tensions a. In the BL case, fitting the n = 1,2,3 data by 

V + U/n + W/n'^ led to a downwards adjustment of V by 1.25 %. The baryonic tensions 

V are listed in Table |2l Motivated by the results of the fits in the mesonic case, the 
baryonic tensions extracted using only n = 1, 2 data were again also endowed with a 
3 % downwards systematic uncertainty. 



^The geometry labels BL, BT, BS originate from the author's subjective picturing of the respective 
configurations as "L-shaped" , "T-shaped" and "symmetric" . 



9 



Geometry 
label 


Static quark 
positions 


Scale factors 
n used 


Baryonic tension in 
lattice units Va? 


BL 


(0,0,0) 
n- (a, 0,0) 
n- (0,a,0) 


1,2,3 


1.518 ±0.002 


BT 


n- (0,a,0) 

n- (a, 0,0) 

n- (-a, 0,0) 


1,2 


9 9Q+0.002 


BS 


n- (a, 0,0) 
n- (0,a,0) 
n- (0,0, a) 


1,2 


o o7fi+0-005 
Z.Zm_o.o73 



Table 2: Baryonic geometries and associated baryonic tensions. 

A law predictions for these baryonic configurations are quite straightforward, since 
the only relevant bond string tensions are o"mioo/2 and auiio/'^ (the reader is reminded 
that the A model assumes half the mesonic string tensions to be associated with the 
bonds). The predictions are 



Vtsa' 



2(aMiooaV2) + V2{aMiwa^ /2) = 1.407 ± 0.001 
2y2(aMiioaV2) + 2(aMiooaV2) = 2.0485 ± 0.0015 
3v/2(aMiioaV2) = 1-924 ± 0.002 



(13) 
(14) 
(15) 



On the other hand, to arrive at Y law predictions, one must perform a potential energy 
minimization as a function of the bond junction position; the interplay of varying bond 
lengths and associated string tensions may favor a different bond junction than the 
one which would be favorable if the string tension were perfectly isotropic. Since the 
string tension is enhanced along diagonal directions, the optimal bond geometry will be 
more aligned with the lattice axes than in the isotropic case. Exploring bond junction 
positions in the spatial 1-2-plane for the BL and BT configurations and minimizing 
the associated potential energy for the range of angular dependences cr(a) depicted in 
Fig. one arrives at the Y law predictions 



V^ra' 



V, 



BL" 
BT"' 



1.523 ±0.008 

9 9e-+0.02 



(16) 
(17) 



Specifically, the expectation is determined by using for a (a) the fit to the measured 
string tension values displayed in Fig. |2l whereas the error bars are determined by using 
the most extreme a{a) dependences allowed for by Fig. |21 

Treatment of the BS geometry is more complicated. If one assumes the 60" discrete 
rotational symmetry of the configuration to be unbroken, then the bond junction position 
is located in the (1, 1, 1) direction viewed from the origin and one can straightforwardly 
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use the range of angular dependences cr{P) allowed for by Fig.|21to analyze this restricted 
case, in complete analogy to the treatment of the BL and BT configurations above. 
Searching for the minimal potential energy under this assumption to arrive at a first 
estimate regarding the Y law prediction for the BS geometry, one arrives at 



vis«' 



= 2.271+^-^^^ . (18) 

(1,1,1) ""-1^ ^ ^ 



On the other hand, there is no guarantee that the bond junction will indeed be located 
in the (1, 1, 1) direction viewed from the origin. It is therefore necessary to explore the 
whole three-dimensional range of possible bond junction positions; however, the associ- 
ated potential energies in general can only be given if the general angular dependence 
of the string tension is supplied. To arrive at an estimate for this dependence, the par- 
ticular cases a{a) and (j(/5) displayed in Fig. El which essentially give the string tension 
at 0° and 45° longitude on a sphere (if one additionally uses (t(/3 = 90°) = cr(a; = 45°)), 
were used to construct an interpolation of the string tension for other longitudes ip. To 
be specific, the form (T{(p) = (1 — sm{2(p))a{ip = 0°) + sm{2ip)a{ip = 45°) was assumed, 
which also reasonably well fits the additional information that the dependence on the 
longitude ip at the "equator", i.e., connecting a{a = 90°) with o"(/5 = 90°), must again 
take the same functional form as a{a). This interpolation was carried out both for the 
fits to the measured string tension values in Fig. |21 as well as for the upper and lower 
extremes, in order to again arrive at an error estimate for the Y law generated. On 
this basis, exploring the whole three-dimensional range of possible bond junction posi- 
tions indeed yields a slightly more favorable asymmetric bond junction position and an 
associated Y law prediction of 

V^sa' = 2.25+°:?^ (19) 

which only deviates little from ()18j) : in analogy to above, the error bars are based on 
minimizing the potential energy using the yj-interpolations of the extreme a{a) and cr{l3) 
dependences allowed for by Fig. |21 Thus, the systematic error associated with restricting 
the range of bond junction positions in arriving at (fTH|) appears to be comparatively 
minor. 

7 Conclusions 

Carefully taking into account the anisotropy in the string tension introduced into the 
random vortex world-surface model by the hypercubic treatment of the vortex surfaces, 
one obtains an unambiguous characterization of the behavior of the baryonic potential. 
Assembling the measured potentials and the corresponding A and Y law predictions 
from the previous section, cf. Table El it is evident that the A law does not agree with 
the measured baryonic potential, whereas the Y law yields rather accurate predictions 
for it. It should be noted that the Y law behavior is seen quite clearly already at 
moderate distances in the vortex model, contrary to full lattice Yang-Mills theory J19j . 
where it is necessary to consider rather large distances. In the full theory, presumably 
the short-distance perturbative A behavior masks the long-distance Y law and pushes its 
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Geometry 


Va' 


V'^a' 


\/^a^ 


BL 


1.518 ±0.002 


1.523 ±0.008 


1.407 ±0.001 


BT 


9 9Q4-U.0U2 
2-2y_o.07 


9 9p:+U.U2 


2.0485 ±0.0015 


BS 


9 27fi+'^-UU5 
^•^'0-0.073 


9 9cr+U.02 
^•^'-'-0.16 


1.924 ±0.002 



Table 3: Comparison of measured baryonic tensions with A and Y law predictions. 

onset to larger separations. In the vortex model, by contrast, such perturbative effects 
are truncated and the Y law induced by the vortices clearly dominates the behavior of 
the baryonic potential. 
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